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Abstract 


We present the complete O(a^) corrections to the Wilson coefficient of the 
unit operator in the zero-recoil sum rule for the B —> D* transition. We 
include both perturbative and power-suppressed nonperturbative effects in a 
manner consistent with the operator product expansion. The impact of these 
corrections on \Vcb\ extracted from semileptonic B —> D* decays near zero 
recoil is discussed. The mixing of the heavy quark kinetic operator with the 
unit operator at the two loop level is obtained. 0{as) corrections to a number 
of power-suppressed operators are calculated. 
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I. INTRODUCTION 


Semileptonic decays of B mesons provide an opportunity to measure the Cabibbo- 
Kobayashi-Maskawa matrix parameter \Vcb\ with minimal theoretical uncertainties (for a 
recent review, see e.g. [|I|). One of the two most popular methods is based on the experi¬ 
mental determination of the zero-recoil B D* transition amplitude by extrapolating the 
experimental decay rate of -B — D* iu to the point of zero recoil momentum, where the 
invariant mass squared of leptons is = {Mb — The hadronic B D* transition 

amplitude for this kinematics is written as: 


(Z)*|c7^756|i?) — 2 Fb*\JMbMb* e* , (1) 

where e is the polarization vector of the D* meson. The zero-recoil form factor Fb* is 
calculable in the short-distance perturbative expansion up to terms i^QCD/mc,bf. These 
nonperturbative corrections cannot be evaluated in a model-independent way at present; 
they are expected to be about —8% |^. 

Existing estimates of the long-distance strong interaction corrections to Fb* are based 
on the sum rules for heavy flavor transitions |^,|^. They relate certain sums of the transition 
probabilities to expectation values of local heavy quark operators in the decaying hadron. 
The zero-recoil sum rule for the spatial components of the axial current can be written in 
the form: 


Fd*\^ + T. I^excr 

e<ii 




'b7r(h') 






ml 





( 2 ) 


The functions and are short-distance (perturbative) coefficient functions, mg 

denotes heavy quark masses mb,c and are 5-meson expectation values of the kinetic 

and chromomagnetic operators, respectively: 

= 4(^) = . (3) 

By Fexc we generically denote transition form factors between B and excited charm states 
with masses Mb + e. They are related to the appropriate structure function of the B meson: 


\Fn*\^ + T. 

e</i 

\Fe^c\‘^ = 

/TT Jo 

[e,q = 0)de , 

(4) 


wt = 

e = Mb — Mb* — qo , 

(5) 


with the invariant hadronic amplitudes h defined as in 0: 

= i J d'^xe T [jl{x) >(0)} (here 
hi,^^{B\f,,\B). ( 6 ) 
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The derivation of such sum rules and their usage for estimating physical form factors is 
explained in detail in |^. 

Because of short-distance perturbative effects the sum over excited states in the l.h.s. 
of the sum rules does not converge at e ~ Aqcd- Instead, for Aqcd < e < mg one has 
Wi{e,q = 0) ~ as{e)€/mQ. This necessitates introducing a cutoff at some energy /r. The 
cutoff serves as a natural normalization point of the effective low-energy operators, most 
appropriate for use with the sum rules. The scale /i satishes the constraint Aqcd /i <C 
As explicitly indicated in Eq. (^, all coefficient functions are also //-dependent. 

The leading term is of primary importance, since it accounts for the short-distance 

perturbative renormalization of the zero-recoil axial current. It is calculable in the pertur¬ 
bative expansion provided // is large enough and belongs to the perturbative domain. If 
/i is too large, it weakens the constraining power of the sum rules, since the heavy quark 
expansion runs in powers of ///mg. 

The 0{as) perturbative corrections to the sum rules were calculated in (see also 
0). Separate pieces of BLM-corrections 0 contributing to sum rules at order O^PqoI) 
were considered in |^-^. The complete BLM resummation of the unit operator coefficient 
function was carried out in P|. The impact of these quasi-one-loop corrections on 

the sum rules proved to be small when one follows the Wilson approach to the Operator 
Product Expansion (OPE) assuming explicit separation of short-distance and long-distance 
contributions. 

More challenging are the genuine, non-BLM 0{al) corrections. Their magnitude is 
crucial for estimating the actual impact of unknown higher-order corrections. Technically 
the most complicated piece corresponding to virtual corrections to heavy quark currents at 
zero recoil was calculated to order 0(al) in [^. In the present paper the remaining second- 


order corrections to the zero-recoil sum rule for spatial components of the axial current 
are computed and the complete two-loop expression for the perturbative coefficient function 
is given. We also derive the two-loop evolution of the kinetic operator. As a byproduct 
of the analysis, we find 0{as) corrections to coefficients of some power suppressed operators 
in the nonrelativistic expansion, and obtain a similar correction to the coefficient of the 
kinetic operator in the sum rule. The 0{a^) correction to the sum rule for the timelike 
component of the vector current {B —>■ D transition) will be given elsewhere. 


II. 0(al) CORRECTIONS TO 

The most efficient way to determine ^a(/^) was suggested in 0 (see also [|^,||). It relies 
on considering the OPE relations of the type of Eq. (j^) in perturbation theory. The structure 
function wa is then given by the weak transitions amplitudes between initial and final states 
consisting of quarks and gluons. Our aim here is to calculate it to order 0{a'^). 


^In reality, this amounts to use /t ~ several units times Aqcd- The existence of such a scale 
in practice is the criterion for applicability of the heavy quark expansion to charm quarks at a 
quantitative level. 
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The l.h.s. has an elastic contribution b ^ c for on-shell quarks and the continuum 
contribution from b ^ c + gluon and 6 —> c -f 2 gluons. The elastic contribution is equal 
to which was calculated in [pT[ . The inelastic part of the structure function 

to order Oi^al) is calculated in the present paper as an expansion in fi/mq. In order to 
determine perturbative coefficients in the sum rules, we need the inelastic part only through 
the leading, second order in /i/mg, because as far as nonperturbative effects emerging from 
local higher-dimension operators are concerned, we account explicitly only for l/mq terms.|^ 
Therefore, in order to calculate with 0{al) accuracy, one has to calculate the 

inelastic part of the structure function the perturbative correction to the coefficient 
function of the kinetic operator and the perturbative expectation value of the kinetic 

operator to the necessary order in as- The expectation value of the chromomagnetic operator 
vanishes in the perturbative expansion to leading order in l/mq. 

Through order the perturbative contribution to the l.h.s. of Eqs. (H), (H) has 

the form (the Feynman diagrams for the inelastic part are shown in Fig. 1): 


T yTy»“'(£,o)d£ = (,,!j 


(2 loop) \ 2 


TT \ VT / 



( 7 ) 


where M is a normalization point for the strong coupling constant in the MS scheme. De¬ 
noting X = mc/mf,, we get: 


Af = + 


= CpAi + CAAi + TuN^At 


Aj{fi) = In (x) 


13 


3 2\ 25 11 25 2 




( 8 ) 

(9) 

( 10 ) 


A»=ln(a;) 

917 949 917 , 


In 


2p 


TUr 


So 111 (M" 

-TT^ H- In — 

24 6 l2/i 


720 1080 


-X 


720^ ’ 


aA/ N 77 89 77 2 2 f 

180 270 180 3 ^ \2^i^ 


( 11 ) 

( 12 ) 


^The l/nig corrections were also calculated Q. In this paper, however, we limit our consideration 
to the leading 1/mg nonperturbative corrections. 
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Although the sum of the elastic b ^ c transition probability and the inelastic excitations 
is more infrared-stable than they are separately, it is still not completely free of power 
suppressed infrared contributions. Its infrared part is given by the expectation values of the 
local operators in the r.h.s. of the sum rule (H). Therefore, by calculating the expectation 
values of the local operators in perturbation theory and accounting for them in Eq.(^, we 
eliminate the contribution of the infrared domain from the Wilson coefficient of the unit 
operator 

The expectation value of the kinetic operator in perturbation theory can be determined 
by using the sum rule for spatial components of the vector current at zero recoil in the heavy 
quark limit 


1 



(e) de 









(13) 


The excitation probability (the perturbative structure function) wY (e) is calculated using 
the same technique as for the axial current. To apply this sum rule, however, we need the 
coefficient functions 'with 0{as) accuracy. For this purpose we perform a nonrelativistic 
expansion of the vector current accounting for the 0{as) corrections. The result of this 
calculation reads: 


C76 |g-=o = {-uiD + V [a X i)]) ifb , 
1 f 1 + X 


(14) 


u = 


mb [ 
1 1 


2x 

- X 


1-CA( ^ 1 ln(x) + 1 


V = 


mb 2x 


l-C, 


TT \4(1 — x) 
as ( 3(1 + x) 


+ C, 


Us ln(x) 
27r 1 — X 


TT \4(1 — x) 


ln(x) -I- 1 


This is the Cl(Q; 5 )-corrected form of the expansion for the vector current given in [Q, 
Eq. (181). The coefficients u and v are obtained by evaluating one-loop graphs shown 
in Fig. 2 in the linear approximation in pb^c/f^Q- Using Eqs. (^) one determines the nor¬ 
malization of the state produced by the vector current: 


^^|(n|c7i6|5)|2 

n 


(n^ -1- 2n^) pI — (n^ — 2uv) Pq 


(15) 


As was shown in Ref. [§], this normalization is nothing but the sum rule of interest. We get 


Stt ^ 




1 — -X -f x" 
o 

1 2 

-X -1- x" 

3 3 


2Cf— 


a, 

TT 

2Cf— 

TT 


1 - -X 4- x" 
o 


31-4X, / 10 

-|- --Inx 1-X -|- X 

X V 9 


41 


1 2 

-X -|- X" 

3 3 


1 ,'12 3 2 

Inx ( - H —X -|- -X 

s. 4 3 4 


(16) 

(17) 


The same expression for can be also obtained by considering the sum rule for a slowly 
moving b quark. Below, this method is applied to derive the Wilson coefficient of the kinetic 
operator which enters the axial sum rule [cf. Eqs. (^OD -(^Tl)]. 
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In order to get the expectation value of the kinetic operator in perturbation theory and, 
therefore, its mixing with the unit operator to 0(ag) accuracy, one has to evaluate the 
l.h.s. of the sum rule (^) to second order in perturbation theory through terms fi^/rriQ. 
Since the chromomagnetic operator does not mix with the unit operator, the perturbative 
contribution as a whole should be identihed with the kinetic operator. 

Performing this calculation, we get, with the 0(0^) accuracy: 



as{2fi) 2 
Up - fi 


+ Cf 


TT 





Ca - ^TrNl 

y 



(18) 



The last term in Eq. (^) represents the non-BLM contribution. Note that the term 
is absent. It means that the hrst-order mixing in the Abelian theory without light flavors 
is not renormalized by effects of higher orders. We note that this fact actually holds to all 
orders in perturbation theory (see Ref. 0). 

Finally, we calculate the coefficient function ^ 7 r(/i) of the kinetic operator in the axial 
sum rule Eq. (^) to order 0{as). The result is as follows: 


ad , 

IT \mQ J 


( 20 ) 


^ (l + ^3) + 


(1) 2 2i d2/i\ 9-17x + 31x 


^ 7x^ , 1 + 22x + 

-In X H— 


18 


( 21 ) 


The last relation is obtained by considering the zero-recoil sum rule perturbatively, in the 
hrst order in ck^ for the initial b quark moving with spatial momentum |p | -C rrib, rric- Let 
us note that for mb 7 ^ me, the zero recoil condition ^ = 0 implies a change in the spatial 
quark velocity. This change leads to gluon bremsstrahlung. The virtual corrections become 
infrared divergent, this divergence, as usual, is compensated by the contribution of the real 
gluon emission. This cancelation, however, brings in an explicit logarithmic dependence of 

on /i. 

We note that this logarithmic dependence is not quite usual. Although the kinetic 
operator has a vanishing anomalous dimension, its coefficient in the sum rule contains the 
logarithm of the cut-off parameter /i at order Because of the power mixing of the kinetic 
operator with the unit operator, the coefficient of the latter has a similar logarithm in the 
power suppressed term We must emphasize, however, that the log/i in Eq. (^T]) 

does not originate from the dependence of the coefficient function on the normalization point 
used for the kinetic operator, but rather from the explicit dependence of the observable we 
consider (axial sum rule) on fi. If we introduce a normalization point v of the operator as 
an independent parameter not equal to p, there will be no log v dependence in and only 
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as log /i will remain. A transparent physical picture underlying this fact will be discussed in 
a separate publication. 

Equations (|^, ([T8| ) and (pOD-(pI|) combined with the known result for allow us 

to obtain to order We note, however, that our perturbative expressions were given in 
terms of the pole masses of the heavy quarks as they appear in this order of perturbation 
theory. To get rid of spurious l/mg infrared effects associated with the pole masses, we 
have to switch to the short-distance masses which have concrete numerical values 

at a given /r, independent of the order of perturbation theory. To order l/rriq this change 
affects only the term 7]\. Unless this is done, the sum rules are formally inconsistent since 
the perturbative expansion of tja has an infrared piece of the order of as{mQ)AQCD/mQ ||TT| . 

In principle, the short-distance masses can be dehned in different ways. Since through¬ 
out this paper the renormalization scheme with the cutoff over the excitation energy is 
implemented, we use the 0{as) relation |]T2 


d»^Q(h) 

d/i 


Then we replace 


\Va 


(2 loop) I^^pole^ 12 


Q I 


-2C], 


= -c, 


^ ^ 2a: 


as 


TT 


^ + — + 0 

,3 tuq 



( 22 ) 


In a: -|- 1 -|- a: 


X 



(1 x) 



(23) 


The hnal result for the coefficient function ^a{a) of fh® operator is obtained by 
combining all terms calculated above. Following |1^, we express the result in terms of 
as even though the rric ^ mb symmetry arguments (which motivate this choice) 

do not apply for ^a{a) because an additional momentum scale is present. 

Collecting all pieces together and neglecting terms we obtain for ^a{a)- 


= (»f “-VqIa*)))' + cAA 


m. 


+ C^F 
= -2 


TT 

2x 




mr 


^A + Cf 


(7)' 


A — X 

42) _ n aF , ri /-A 


In a: -|- 1 -|- a: 

CF = CfC + CaH + T„NF’i, 


Ct = xll 


X 


I'ln 

\mcJ 


a;(17 + 5a; + 4a;2) , 25 8 25 o 

-;---In X- X - X , 

6(1-x) 18 9 18 ’ 


(24) 

(25) 

(26) 

(27) 
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y ^ 18 ^ 


, 203 1859 

Inx H-1- 

80 1080 
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' 2x 
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In 
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203 

W 

77x 


-X 


71x2 

yr 


(28) 

(29) 
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The function to order can be found in 

Since the BLM part of the corrections was discussed in the literature in detail [0, we 
single out the genuine non-BLM part which is defined as the value of the full second-order 
coefficient 02 at ■ 


Q. / Q; \ ^ 

Uif^) = 1 + ai rngifi)) — + 02 (/X, rngin)) — + ... 

71 \ 71 / 


a. 


TT 


= 1 + ai in, mgifi)) — + (/x,mQ(/i))/?o + 02 (/^)’^q(/^)) 


a< 


71 


+ 


11 4 

f^O = ~ 


(30) 

(31) 


Note that 0 ° does not depend on the convention for the normalization point of the strong 
coupling. Its value is shown in Fig. 3 as a function of /x/roc for three values of rric/mh = 
0.2, 0.25 and 0.3. We denote the corresponding non-BLM second-order shift in ^yi(/i) by 
A°(/x): 


A°(/i) = alifi, mg{n)) . 


III. \Vcb\ determination and zero-recoil sum rules 


Let us now turn to the application of our results for the determination of |14fe|- First, we 
note that the net impact of the non-BLM a"^ corrections on the sum rules is rather small. 
Taking a reasonable value /x/ruc = 0.5 and rric/mb = 0.25 we get ~ —0.7. Assuming 
as = 0.22 to 0.27, the absolute shift A°(/x) constitutes —0.003 to —0.005 which translates 
into —0.0015 to —0.0025 decrease of the short-distance perturbative renormalization of the 
zero-recoil axial current. For any reasonable choice of the parameter /x in the sum rules 
consistent with using the l/rric expansion, this effect is well below 1%. It should be noted 
that since in the Wilson OPE the infrared domain is completely excluded from the coefficient 
functions, the effective coupling cannot become large. 

The perturbative correction to the coefficient of the kinetic operator (Eq. ([2l|)) is strongly 
suppressed. The actual value of this coefficient changes only by several percent. We did 
not calculate the corresponding effect in the chromomagnetic term; due to the anomalous 
dimension of this operator it depends on the specific renormalization procedure. In view 
of the result for the vector sum rule (see Eq. ([T7|)) we do not expect this correction to be 
significant either. 

The axial sum rule (j^) allows one to get a QCD-based estimate of the combined effect 
of the perturbative and nonperturbative effects on the B —> D* zero-recoil form factor: 


Ed* I — 


^ ^ ^7r(/x) 2 / ^ Cg(/^) 

mi mi 


11/2 


/^G (/^) ^2 I 




+ o 




m 


(32) 


Q, 


It is seen that plays the role of a short-distance renormalization factor of the weak 

current, while the remaining terms in square brackets yield 1/m^ long-distance corrections 












to the form factor. The quantum-mechanical meaning of each of the three power terms, 
conhrming this formal conclusion, was elucidated in [Q (see also |]l|). 

With the good perturbative control over the short-distance part of Fd*, the main uncer¬ 
tainty dominating theoretical predictions for the form factor resides in power corrections. 

In Ref. p] higher order BLM-type corrections to the sum rule were analyzed. It was 
shown that, assuming that the running of the QCD coupling as below the charm mass is 
a valid practical concept, it is necessary to perform a resummation of the leading BLM 
corrections in order to arrive at a meaningful numerical result. Within the Wilson approach 
to OPE, the overall impact of the BLM corrections appears to be moderate. The typical 
BLM-resummed value for at p/mc — 0.5 appears to be near 0.99. 

On the practical side, excluding the infrared domain notably improves the accuracy 
of the perturbative calculations for the charm quark, at least in the context of the BLM 
calculus. We recall that the perturbative zero-recoil factor rjj^ has an intrinsic uncertainty ~ 
{Aqcd/t^cY due to infrared renormalons. By calculating and taking into account 

fA^/rriQ terms, the (AgcD/’^c)^ uncertainty is eliminated. As long as we do not account 
for l/rriQ terms explicitly, the perturbative expressions still have an infrared renormalon 
uncertainty of the order of (Aqcd/’^c)^- Using the same overall normalization for both 
cases one has 



Performing a simple estimate, one hnds that the size of the 1/m^ uncertainty in r]A becomes 
quite signihcant for Aqcd — 250 MeV. This shows up in the BLM corrections which become 
quite large already at the lowest orders. Shifting the uncertainty down to (Aqcd/'"2-c)^, 
we significantly reduce it. 

The analysis of the 0(ag) corrections presented in this paper shows that the genuine two- 
loop effects are quite small. Therefore, our numerical conclusions do not differ in practice 
from the estimate of the zero-recoil axial form factor given in ||l|] 


D* 


~ 0.91 - 0.013 


fj.1 - 0.5 GeV" 


± 0.02excit ± O.Olpert ± 0.025i/m3 


(34) 


0.1 GeV^ 

For further improvements, one has to bring in a new dynamical input yielding the magnitude 
of long-distance l/m^ and 1/m^ corrections more precisely. 

Experimental data on the small-recoil B —>■ D* decay rate are not fully conclusive yet. 
Nevertheless, the value of iWftI extracted from the exclusive B D* transitions using Fd* ^ 
0.9 seems to be in a reasonable agreement — already within experimental uncertainties — 
with the results obtained from inclusive semileptonic decay width Psi(il —>■ Xclvi). The 
estimate of the complete 0(a'^) correction in rsi(5 —>■ XJui), recently presented in Ref. |T4 


is another example of the theoretical progress in the perturbative treatment of one of the 
most important heavy quark decays. 


IV. CONCLUSIONS 

We have calculated the complete (P(a^) corrections to the zero-recoil heavy quark axial 
sum rule which is used to evaluate the zero-recoil B D* form factor. The calculations 
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performed here are a necessary supplement to the results of Ref. The use of the sum 


rules allows one to incorporate both O^a"^) and power-suppressed nonperturbative effects in 
a consistent way. The genuine (non-BLM) two-loop corrections are shown to be relatively 
small and under good theoretical control. We note that the uncertainty associated with 
these effects had not been reliably estimated previously. This theoretical uncertainty is now 
eliminated. 

As an important theoretical conclusion, we emphasize that a consistent implementation 
of the Wilson OPE separating short- and long-distance contributions, is feasible even when 
highly nontrivial complete 0(al) corrections are taken into consideration. In principle, 
this procedure does not bring in additional complications, compared to purely perturbative 
calculations performed without an infrared cutoff. The Wilson approach, on the other 
hand, allows one to operate with well dehned notions of short-distance and long-distance 
effects. From the practical viewpoint, it allows one to decrease signihcantly the uncertainty 
in perturbative coefficients in beauty decays. This essential reduction originates mainly from 
the BLM-type corrections; the genuine two-loop effects are less radically changed. 

We note also, that only in the framework of the Wilson OPE approach to QCD is 
it possible to preserve a number of exact inequalities formulated for hadrons containing 
a heavy quark. They exist in a simplihed quantum mechanical treatment which ignores 
the peculiarities of the held-theoretical description. These inequalities are important in 
constraining those parameters of the heavy quark expansion which are not yet measured in 
experiment. 

As a byproduct of our analysis, we obtained the complete two-loop perturbative evolution 
of the kinetic operator, Eq. (0fl 

We also calculated 0{as) corrections to the coefficient function of the kinetic operator in 
the axial sum rule. The nonrelativistic expansion of c^b at zero recoil and the corresponding 
sum rule were obtained with 0{as) accuracy. Numerically, we found an overall short-distance 
renormalization of the zero-recoil current to be very small, close to the estimates of [^,^,|^ 
and rather different from the value used in other analyses of Fd* where long-distance 1/m^ 
corrections were addressed. 


ACKNOWLEDGMENTS 

N.U. gratefully acknowledges discussions of the OPE subtleties with the TPI fellows 
M. Shifman, M. Voloshin and A. Vainshtein, and valuable perturbative insights from 
Yu. Dokshitzer. He also thanks I. Bigi for encouraging interest and R. Dikeman for com¬ 
ments. This work was supported in part by DOE under the grant number DE-FG02- 
94ER40823, by BMBF under grant number BMBF-057KA92P, and by Graduiertenkolleg 
“Teilchenphysik” at the University of Karlsruhe. 


^The BLM part was analyzed to all orders in Q. 


10 




REFERENCES 


[1] I. Bigi, M. Shifman, N.G. Uraltsev, Preprint TPI-MINN-97/02-T, March 1997, hep- 
ph/9703290, to be published in Ann. Rev. Nucl. Part. Sci. 

[2] M. Shifman, N.G. Uraltsev and A. Vainshtein, Phys. Rev. D51, 2217 (1995), erratum: 
ibid. D52, 3149 (1995). 

[3] I. Bigi, M. Shifman, N.G. Uraltsev and A. Vainshtein, Phys. Rev. D52, 196 (1995). 

[4] J. G. Korner, K. Melnikov, O. I. Yakovlev, Zeit. f. Physik C69, 437 (1996). 

[5] S. Brodsky, G. Lepage and P. Mackenzie, Phys. Rev. D28, 228 (1983). 

[6] M. Neubert, Phys. Lett. B341, 367 (1995). 

[7] P. Ball, M. Beneke and V. Braun, Phys. Rev. D52, 3929 (1995). 

[8] A. Kapustin et al, Phys. Lett. B375, 327 (1996). 

[9] N.G. Uraltsev, Nucl. Phys. B491, 303 (1997). 

[10] A. Gzarnecki, Phys. Rev. Lett. 76, 4124 (1996); 

A. Gzarnecki and K. Melnikov, hep-ph/9703277, submitted to Nucl. Phys. B. 

[11] N. G. Uraltsev, Int. J. Mod. Phys. All, 515 (1996). 

[12] 1. Bigi, M. Shifman, N. Uraltsev and A. Vainshtein, Preprint GERN-TH/96-191, hep- 
ph/9704245, submitted to Phys. Rev. D. 

[13] N. G. Uraltsev, Mod. Phys. Lett. AlO, 1803 (1995). 

[14] A. Gzarnecki and K. Melnikov, Phys. Rev. Lett. 78, 3630 (1997); 

A. Gzarnecki and K. Melnikov, hep-ph/9706227, submitted to Phys. Rev. D. 


11 



FIGURES 



FIG. 1. Examples of inelastic contributions to the zero recoil sum rules. 
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FIG. 2. One-loop diagrams determining the a^-corrections to the nonrelativistic expansion of 
the vector current in the small velocity kinematics. 
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